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We show that non-linear optical structures involving a balanced gain-loss profile, can act as unidi- 
rectional optical valves. This is made possible by exploiting the interplay between the fundamental 
symmetries of parity ("P) and time (T), with optical nonhnear effects. This novel unidirectional 
dynamics is specifically demonstrated for the case of an integrable "PT-symmetric nonlinear system. 
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I. INTRODUCTION 



Transport phenomena and in particular directed trans- 
port are at the heart of many fundamental problems in 
physics, chemistry and biology [jLj. At the same time 
they are also of great relevance to technological appli- 
cations based on a variety of transport-based devices 
such as rectifiers, pumps, particle separators, molecu- 
lar switches, and electronic diodes and transistors. Of 
special interest is the realization of novel classes of in- 
tegrated photonic devices that allow one-directional flow 
of information-e.g. optical isolators [2|. Currently, such 
unidirectional elements rely mainly on the Faraday ef- 
fect, where external magnetic fields are used to break 
the space-time symmetry. This in general requires ma- 
terials with appreciable Verdet constants-typically non 
compatible with light-emitting wafers [5]. To anticipate 
these problems, alternative proposals for the creation of 
optical diodes and isolators have been suggested in re- 
cent years. Some representative examples include the 
creation of optical diodes based on asymmetric nonlinear 
absorption [3], second harmonic generation in asymmet- 
ric waveguides p], nonlinear photonic crystals [5], and 
photonic quasi-crystals and molecules [6 . 

In this paper, we propose a new mechanism for uni- 
directional optical transport based on configurations in- 
volving nonlinear optical materials with Parity (V) and 
Time (T)-reflection. This is possible by judiciously in- 
terleaving gain and loss regions, in such a way that the 
(complex) refractive index n{x) = nji{x) +ijni(x) pro- 
file satisfies the condition n*{—x) = n{x). A first ex- 
perimental realization of such (linear) arrangements has 
been recently reported in Refs. |71 ^ where a VT dual 
coupled structure was fabricated and the beam dynam- 
ics was investigated. Here we show that the interplay of 
non- reciprocal dynamics arising from T-'T-symmetry 8J, 
and self-trapping phenomena associated with Kerr non- 
linearities [51 , can mold the flow of light in a surprising 
way. Such novel directed dynamics can be exploited in 
the realization of a new generation of optical isolators or 
diodes. 



Even though the validity of our arguments can 
be demonstrated for a variety of non-linear VT- 
configurations, below, we will highlight its basic princi- 
ples, using the simplest possible arrangement, consist- 
ing of two PT-coupled waveguide elements with Kerr 
nonlinearity of strength x- Each of the waveguides is 
single-moded- one providing gain and the other an equal 
amount of loss (see Fig. [2]). We have obtained the phase 
diagram in the x ~ 1 plane for which our system acts 
as an optical diode, and we have identified the minimum 
propagation length needed, in order to achieve this uni- 
directional functionality. Detail numerical simulations 
support our theoretical predictions. 

This paper is structured as follows. In Sec. II an 
overview of the linear T'T-symmetric dimer is presented. 
The nonlinear T'T-symmetric dimer will be introduced in 
Sec. Ill, where the equations of motion are given in terms 
of Stokes parameters. In subsection III. A, we present 
both our theoretical and numerical results on the dynam- 
ics of the non-linear T^T-symmetric dimer. In subsection 
III.B we calculate the critical value of the non-linearity 
for which diode action is possible. Finally we will draw 
our conclusions in Sec. IV. 



II. LINEAR PT-SYMMETRIC DIMER: 
OVERVIEW 



AN 



In this section we will briefly review the basic prop- 
erties of the linear T'T-symmetric dimer [71 |H1 E]. In 
integrated optics this simple VT element can be realized 
in the form of a coupled system, with only one of the two 
parallel channels being optically pumped to provide gain 
7 for the guided light, whereas the neighbour arm expe- 
riences equal amount of loss (see Fig. 1). Under these 
conditions, and by using the coupled-mode approach, the 
optical-field dynamics in the two coupled waveguides are 
described by the following set of equations: 



dz 
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ili^i = 0; 
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where ^^1,2 are modal electric field amplitudes in the am- 
plifying (Eq.[I^) and lossy (Eq. ^) waveguide channels, 
z represents a dimensionless propagation distance- nor- 
malized in units of coupling lengths, and 7 is a scaled 
gain/loss coefficient, also normalized to the coupling 
strength. 

The Hamiltonian corresponding to the linear problem 
of Eq. ([I]), is written as: 



H = 



17 
-1 



-1 



(2) 



and commutes with the combined VT operator. A sur- 
prising result associated with this class of problems is the 
possibility that such a VT symmetric Hamiltonian H can 
have an entirely real energy spectrum, despite the fact 
that, in general, they are non-Hermitian [71l51 [TT| - [2n] . For 
the specific example of the non-hermitian Hamiltonian of 
Eq. (I2]), a direct diagonalization gives the following set 
of eigenvalues: 



A± = ±yr^ 



Y 



(3) 



which are real as long as the gain/loss parameter 7 is 
smaller than some critical value, j-pT — 1 {exact VT- 
symmetric phase). As the gain/loss parameter 7 in- 
creases above jvT: the eigenvalues becomes complex 
(broken T'T-symmetric phase). The corresponding eigen- 
vectors of the Hamiltonian Eq. ^ are 
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In the exact T^T-symmctric phase, both the H and 
VT operators share the same set of eigenvectors. In this 
regime, the mode intensity is symmetric with respect to 
the mirror axis of the two waveguides. As 7 increases 
above j-pT the eigenfunctions of H cease to be eigenfunc- 
tions of the T^T-operator, despite the fact that H and the 
PT-operator still commute. This happens because the 
■pT-operator is anti-linear, and thus the eigenstates of H 
may or may not be eigenstates of VT- In the broken VT- 
symmetric phase, the spacial distribution of the modes is 
asymmetric, one of them living predominantly in the am- 
plifying waveguide and the other in the lossy one. At the 
phase-transition point 7 = j-pT the two eigenfunctions 
and their corresponding eigenvalues coalesce leading to 
an "exceptional" point singularity pT] , 

The beam dynamics associated with Eq. (fTl) were in- 
vestigated theoretically in Refs. [3 [TTl [T^] while direct 
measurements were performed in [7l |8]. These authors 
recognized that as the gain/loss parameter 7 reaches 
7-P7-, the total beam power starts growing exponentially, 
while for 7 < 7^7- power oscillations are observed (see 
Figs. lb,c). The most dramatic effect in the beam evo- 
lution is the appearance of non-reciprocal wave propaga- 
tion (see Figs. Ic-f). Specifically, the beam propagation 
pattern differs depending on whether the initial excita- 
tion is on the left or right waveguide. This is contrasted 
with the 7 = case (Fig. la,b), where the beam propa- 
gation is insensitive to the initial condition. 




FIG. 1. (color online) Beam propagation in two coupled 
linear waveguides. For the parameters of the simulation (we 
use normalized coupling units) , the spontaneous "PT-breaking 
take place at 'yvT = 1- In all cases, left (right) panels cor- 
respond to an initial excitation at the left (right) channel. 
(a,b) A passive system corresponding to 7 = 0. The propa- 
gation is reciprocal; (c,d) 7 = OA'y-pT corresponding to the 
exact "PT-phase. A non-reciprocal beam propagation is ev- 
ident. Although the dynamics is non-Hermitian, the evolu- 
tion is "pseudo-unitary" and the total beam power remains 
bounded; (e,f) For 7 — 1.67737-, beam power grows exponen- 
tially (vertical scale is logarithmic) in both waveguides, while 
the beam propagation is again non-reciprocal with respect to 
the mirror axis of the two waveguides. Waveguides are color- 
coded, indicating balanced gain (red) and loss (green) regions. 
Grey-color waveguides indicate a passive (7 = 0) system. 



III. NON-LINEAR PT- SYMMETRIC DIMER 

We begin our analysis by providing the mathemati- 
cal model that describes optical wave propagation in a 
Kerr nonlinear VT symmetric coupled dual waveguide 
arrangement (see Fig. [2]). The two modal field ampli- 
tudes are governed by the evolution equations: 



dz 



+ i>2-il'^i+x\^i?'^i = 0; (a) 
+ V^i+*7V'2+xlV'2pV'2 = 0; (6) 



(5) 



where x is the strength of the Kerr-nonlinearity. 

Equations (Im) can be rewritten in terms of the (real) 
Stokes parameters Si — ^p^&iip, where ai{i = 0, 1, 2, 3) de- 
note the Pauli spin matrices j22j . In this representation, 
the total field intensity is given by 5o = If^iP + |V'2p, 
S3 — IV'iP ~ l^2p is the intensity imbalance between 
the two waveguides, while 5*1 = V'i'02 + '0i'02 and S2 = 



form: 



"01 "02)- In this representation Eqs. (Is]) take the 
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where we have introduced the two real vectors E = 
(0,0,27) and B = (2,G,xS'3), and the 3-dimensional 
Stokes vector S = (5*1, ^2, 6*3). We note that the con- 
dition Sq — S ■ S — is always satisfied. It is worth 
mentioning that Eqs. m, are identical to the equation 
of motion of a relativistic negatively -charged particle 
with zero mass, in a pseudo-electromagnetic field {E , 13), 

where {So,o) represents the energy and 3— dimensional 
momentum of the particle, while the propagation dis- 
tance z has the role of the time. 

Nonlinear PT-symmetric optical coupled systems can 
be realistically synthesized on semiconductor wafers- 
known for their high Kerr-like nonlinearities 10 ■ As in 
Ref. [8 , coupling lengths as low as L^ = 1mm can be 
obtained, in which case a gain/loss level below zL30cm~^ 
(readily available in such materials) will suffice to keep 
the arrangement in the VT phase. In addition, critical 
switching (x ~ 1) can also occur at milliwatt power levels 
in multi-quantum well configurations. 

A. Dynamics 




FIG. 2. (color online) Beam propagation in two coupled non- 
linear waveguides with non-linearity strength x and a complex 
7'T-symmetric refractive index profile. Waveguides are color- 
coded, indicating balanced gain (red) and loss (green) regions 
(7 = 0.1). Left columns correspond to an initial excitation at 
the gain waveguide port, while right columns correspond to 
an initial excitation at the lossy waveguide: (a,b) The non- 
linearity X ~ 1-9 is below the critical value Xd ~ 3.37 while 
for (c,d) the non-linearity strength x = 8 is above. 

For 7 = 0, Eqs. ^ admits two constants of mo- 
tion: the total input power Sq and the total energy 



H — (x/2)S'| -I- 251. These two constants allow for an 
analytic solution of the Stokes vector S in terms of el- 
liptic functions [9 . Depending on the initial preparation 
and strength of nonlinearity x, we observe two distinct 
dynamical behaviors. For example, if the initial beam of 
total input power 6*0(0) = 1, is prepared in one of the two 
waveguides (i.e. 63(0) ~ ±1), we observe either Rabi os- 
cillations, or self-trapping dynamics [9]. The former case 
corresponds to x < 4 and results in beam oscillations 
between the two waveguides, while the latter case, oc- 
curs for X > 4 and leads to localization of the field (for 
all times) at the waveguide that was initially placed. In 
both cases, symmetric initial preparation will result to a 
dynamics which is reciprocal with respect to the axis of 
symmetry of the two coupled waveguides. 

For 7 7^ 0, neither the energy "H nor the beam power 
So are anymore conserved quantities. Nevertheless, VT- 
symmetry enforces two other constants of motion C, J: 



J 



So 
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xSi-2 
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thus indicating that the system of Eqs. (Is]) is fully in- 
tegrable. Below we will consider the case where initially 
5o(0) = 1, 5*3(0) = ±1, while 5i(0) = 52(0) = 0. In this 
case, the constants of motion, as defined in Eqs. (I?]), take 
the values C± = ±2 and Jj- = 1 =p 77r/x. 

Using C and J, in this particular case we can express 
the components of the Stokes vector in terms of So{z) in 
the following way 



X^i = 2 (1 - cos(xi^)) , (a) 

XS2 = 2sin(x^^), (b) 

XS3 - ±^{xSo{zW - (4sin (^(1 - So{z)))y. (c) 

Substituting the expression for ^3 from Eq. (|8p), to the 
first of the Eqs. (pi), we get that 



± 



Soiz) 



dSo 



27 



So = l 



^{xSoy-{'isin[^{l-So))y 



z. (9) 



Even though the problem is soluble by quadratures, the 
integral in Eq. ^ cannot be evaluated further and thus a 
closed expression for So{z) is not possible. It is therefore 
instructive at this point to gain insight on the proper- 
ties of the dynamics of this VT non-linear coupler by 
numerically solving Eqs. (5]6|. The accuracy of the nu- 
merical integration was checked via the conservation laws 
Eq. ([7]), which were satisfied up to 10"^". 

Examples of the resulting beam dynamics for 7 = 0.1 
and two representative non-linearity strengths x — 1-9 
and X = 8 are reported in Fig. [2^,b and [2]:,d respec- 
tively. In contrast to the 7 = case 9,, now the dy- 
namics is non-reciprocal with respect to the axis of sym- 
metry of the system. While this is true for both values 



^r 




FIG. 3. (color online) Main figure: A semi-logarithmic plot 
of Xd vs. 7. For the immerical evaluation of Xd we have 
integrated Eq. ([5| Inset: The efficiency factor Q vs. nonlin- 
earity strength x, for a fixed gain/loss parameter 7 = 0.1 and 
three different waveguide lengths z — 10, 20 and z = 30. For 
non-linearity strength X ~ Xd ~ 3.4 the isolator reaches its 
optimal efficiency. 



of non-linearity strength x, it is much more pronounced 
for the case of Figs. [2];,d. In this latter case, the output 
field always leaves the sample from the waveguide with 
gain (red-colored) irrespective of the preparation of the 
input beam. At the same time the output beam inten- 
sity at the lossy waveguide approaches zero for longer 
waveguides. It is important to stress that in the case 
of the linear PT-dimer (see Figs le-f) the beam inten- 
sity at the lossy waveguide never goes to zero. Instead 
it increases exponentially, albeit with a smaller prefac- 
tor with respect to the one of the gain waveguide. This 
novel unidirectional propagation of the T^T-symmetric 
non-linear dimer is the key mechanism for establishing 
optical isolators (diodes). It has to be contrasted with 
the corresponding cases shown in Figs. [2k,b where the 
output beam depends on the input state, i.e. an initial 
excitation at the gain waveguide results in an output field 
at the lossy guide and vise versa. 

To quantify the ability of our set-up to act as an opti- 
cal non-reciprocal device, we have defined the efficiency 
factor Q of unidirectional propagation as 



Q(z) = l-|r+,+ (z)-T_,+ (z) 



(10) 



where T±^+{z) = \'ipi{z)\'^ / So{z) is the normalized trans- 
mission coefficient associated with the gain (-I-) waveg- 
uide of length z. In our definition we have always 
assumed that the initial input beam has total power 
So{z = 0) ~ 1, while the beam is launched either in 
the gain (-I-) or in the loss (— ) waveguide. The efficiency 
factor takes values from < Q < 1: a perfect diode cor- 
responds to Q = 1 (since the term inside the absolute 
value in Eq. (10 1 will be zero), while the opposite limit 
of Q = indicates total revival of the field. In the inset of 
Fig. [3] we report our numerical findings for the efficiency 
factor Q as a function of the non-linearity strength x for 



three different waveguide lengths z — 10, 20, and z = 30, 
and for a fixed value of the gain/ loss parameter 7 = 0.1. 
It is clear that an optimal diode is achieved once the non- 
linearity strength x is larger than a critical value Xd- 
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FIG. 4. (color online) The numerically extracted first passage 
distance Zfpd versus the gain/loss parameter 7. The initial 
conditions are chosen to be 5*0(0) = 1 and 5*3(0) — —1. An 
inverse power law is observed. Left inset: The proportionality 
coefficient /(x) is plotted versus the nonlinearity strength x 
for X > Xd- The red line correspond to the best linear fit. 



B. Critical Non-linearity 

Next we present a heuristic argument that aims to esti- 
mate the critical non-linearity strength Xd (as a function 
of 7), above which the VT symmetric non-linear dimer 
acts as an optical diode of high efficiency factor Q — 1. 
To this end we focus our analysis on the temporal behav- 
ior of the total power So{z). In the case of (Rabi-like) 
oscillations So{z) is bounded between a minimum and a 
maximum value. Instead, in the regime where the cou- 
pled system acts as an optical diode, So{z) is bounded 
only from below, while asymptotically it grows in an ex- 
ponential fashion [53] . Using the first of Eq. ^ together 
with Eq. ([8p), and requesting the extrema condition 
dSoiz)/dz = (which is equivalent to 5*3(2:) — 0) to- 
gether with the condition (PSo/dz'^ < for the existence 
of a global maxima, we find that Sq{z) shows oscillatory 
behavior (i.e. Rabi-like oscillations) if the non-linearity 
X is smaller than Xdj given by 



Xd 



27r7 



(11) 



In the main panel of Fig. p] we compare Eq. ( |11[ ) with the 
numerical values found for Xd- The latter has been eval- 
uated via a direct integration of Eq. (|5| for systems sizes 
up to z = 10^. The critical nonlinearity Xd was evaluated 
up to a fourth digit accuracy as the nonlinearity strength 
for which the total power So{z), is bounded. In all cases 
the accuracy of the integration scheme has been guaran- 
teed by requesting that the constants of motion Eq. (It]) 



are conserved with accuracy up to 10~^. A nice agree- 
ment between the theoretical and numerical value of Xd 
is evident for small values of the gain/loss parameter 7, 
while deviations from the theoretical prediction start to 
be visible as 7 approaches the VT transition point (i.e. 
7 = 1) of the linear system. 




FIG. 5. (color online) Dynamics of the rescaled Stokes vari- 
ables F for X = 9 and various gain/loss parameters 7: dashed 
blue line correspond to 7 = 0.157; pink line (*) correspond 
to 7 = 0.15; solid yellow line correspond to 7 = 0.12; and 
red line (■) to 7 = 0.174. The green line (•) corresponds to 
the passive system 7 = with critical non-linearity x = 4, 
where the motion of the trajectory is on the separatrix. The 
trajectory associated with 7 = 0.157 (see red arrow in Fig. B| 
is typical to the cases where Zfpd diverge and correspond to 
the closest one to the separatrix of the passive 7 = system. 

Finally, we investigate the minimal waveguide length 
Zd which is required in order to have a high— Q diode. 
From Figs. [2j;,d we see that the beam evolution follow 
two distinct scenaria depending on the initial conditions: 
if the beam is launched initially at the gain waveguide, 
the propagation is mainly along this channel. If on the 
other hand the beam excites the lossy waveguide, there is 
a minimum propagation distance Zd which is required be- 
fore the light intensity is concentrated in the gain waveg- 
uide. We have found that Zd is proportional to the "first 
passage distance" Zfpd associated with the point that S'3 
becomes zero for the first time. In Fig. [4] we report 
the results of our simulations for Zd ~ Zfpd for various 
x(> Xd) values or input power levels. 

An intriguing feature of Zfpd is the existence of singu- 



larities (peaks in the Zfpd) for some characteristic values 
of the gain/loss parameter 7. To understand the origin of 
these singularities, we have plotted the evolution of the 
Stokes vector S , by making use of the rescaled variables 
F = S /Sq. In this representation, the magnitude \T\ 
remains constant, and thus we can visualize the evolu- 
tion on the Bloch sphere (see Fig. [5]). It should be em- 
phasized that the Bloch trajectories can in general show 
self-intersections, as they are a projection from a higher 
dimensional phase space. One must also distinguish be- 
tween closed orbits and those approaching an asymptotic 
state, as this is in general connected to broken and un- 
broken T'T-symmetry [T3]. We note that closely related 
Bloch dynamics appear in different physical model sys- 
tems like the ones reported in [21j. Our analysis, indi- 
cated that the singularities in Zfpd are associated with 
trajectories that during their evolution, they stay close 
to the separatrix associated to the critical value x = 4 
(transition between Rabi-oscillations and self-trapping) 
of the passive system. 

Leaving aside the issue of the singularities, we have 
found that for all x- values larger than Xd, the first pas- 
sage distance Zfpd follows an inverse power law i.e. 



^^fpd = /(x)/7 



(12) 



where the proportionality factor /(x) is x-dependent. A 
best least square fit allow us to extract the various /(x) 
which is in this case /(x) — —0.6 + 0.51n(x) (see inset of 
Fig. |4]). 



IV. CONCLUSIONS 

In conclusion, we have proposed a new mechanism for 
directed transport in nonlinear optical coupled systems, 
that relies at the interplay between nonlinearity and VT- 
reflection symmetries. More specifically, we have ob- 
served that above a critical non-linearity strength, the 
beam evolution is unidirectional, i.e. the output beam 
remains in the gain channel, irrespective of initial condi- 
tions. Such behavior implies that these systems can be 
used to realize new classes of optical diodes and other 
unidirectional photonic elements. Of great interest will 
be to extend these notions to more involved arrangements 
like nonlinear VT lattices where nonlinear excitations are 
expected lead to even more intriguing phenomena. 
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